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ABSTRACT 

The redshift dependence of the abundance of galaxy clusters is very sensitive to the statistical 
properties of primordial density perturbations. It can thus be used to probe small deviations 
from Gaussian initial conditions. Such deviations constitute a very important signature of 
many inflationary scenarios, and are thus expected to provide crucial information on physical 
processes which took place in the very early Universe. 

We have determined the biases which may be introduced in the estimation of cosmolog- 
ical parameters by wrongly assuming the absence of primordial non-Gaussianities. Although 
we find that the estimation of the present-day dark energy density using cluster counts is not 
very sensitive to the non-Gaussian properties of the density field, we show that the biases can 
be considerably larger in the estimation of the dark energy equation of state parameter w and 
of the amplitude of the primordial density perturbations. 

Our results suggest that a significant level of non-Gaussianity at cluster scales may be 
able to reconcile the constraint on the amplitude of the primordial perturbations obtained 
using galaxy cluster number counts from the Planck Sunyaev-Zeldovich Catalog with that ob- 
tained from the primary Cosmic Microwave Background anisotropics measured by the Planck 
satellite. 
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1 INTRODUCTION 

Inflation is widely accepted as the most elegant mechanism to gen- 
erate the density fluctuations which seeded the Large-Scale Struc- 
ture (LSS) of the Universe we observe today. The simplest standard, 
single field and slow-roll, inflation model predicts, among other 
things, that the primordial density perturbations were nearly Gaus- 
sian distributed (see e.g. |Creminelli|2003[ |Maldacena| 2003 1 |Lyth 



& Rodriguez]|2005| [See^ & Lidsey||2005| [Sefusatti & Komatsu 



2007 1. While such predictions seem to be in good agreement with 
current observations of the cosmic microwave background (CMB) 
anisotropies (e.g. |Slosar et al.|2008^ and large-scale structure (e.g. 
[Komatsu et al.|2011 1, a significant, potentially observable level of 
non-Gaussianity may be produced in many inflationary models. 
Detecting and constraining primordial non-Gaussianities has be- 
come a crucial task in current cosmological studies, since a pos- 
itive detection of non-Gaussianity would rule out a considerable 
number of inflationary models, opening an entirely new window 
into the very early Universe. 

A wide range of observables have been used to probe primor- 
dial non-Gaussianities. While the three-point statistics of the tem- 
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perature anisotropies in the CMB is the most common tool, the 
statistical properties of the large-scale structure, namely the bispec- 
trum and/or trispectrum of the galaxy distribution (e.g.|Sefusatt i &| 
|Komatsu|2()07l|Matarrese & Verdel200 8 Giannantonio et al. IQV^ 
and weak-lensing observations (e.g.'Schafer et al. 2012. |Hilbert et] 
|al.|2012> , as well as CMB-LSS (Tashiro & Ho 2012]l and CMB- 
21cm ( [Takeuchi et al.|20T2 ) cross-correlations have been used to 
the same effect. The evolution with time of the abundance of mas- 
sive collapsed objects such as galaxy clusters (see e.g. [Matarrese] 
|et al.|2000{ [Robinson & Baker|200^ also holds key information 
that could be used to probe primordial non-Gaussianities. In fact, 
we have shown (Trindade e t al.|20T2y that assuming the absence of 
primordial non-Gaussianities, may lead to an apparent discontinu- 
ity in the evolution of the estimated effective equation of state pa- 
rameter with redshift using galaxy cluster counts. In the same spirit, 
we will estimate and quantify the biases which might be introduced 
in the determination of several of the most important cosmological 
parameters using the evolution with time of the galaxy cluster abun- 
dance, if it was wrongly assumed that the initial conditions of the 
primordial density field were Gaussian distributed. We focus our 
attention on the present-day dark energy density, Q^^, the dark en- 
ergy equation of state parameter w (here assumed to be a constant), 
and the present-day root mean square mass perturbations, erg, at 
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the standard 8 /?"' Mpc scale. Primordial non-Gaussianity is usu- 
ally parametrized by a non-linear parameter, /^vl^ which relates the 
bispectrum and power spectrum of the primordial curvature pertur- 
bation jLo Verde et al.|2008[ [Trindade et al.|2012[ l. The parameter 
/wL can be calculated using different prescriptions, each associated 
with a different physical mechanism for the generation of primor- 
dial non-Gaussianities. We will mostly consider the so-called Local 
configuration (see [Trindade et al.|2012| for more details), although 
we will also mention at some point the Equilateral configuration. 



2 GALAXY CLUSTER ABUNDANCE 

The abundance of galaxy clusters, as a function of mass and red- 
shift, contains information which can be used to compare cosmo- 
logical models or constrain the value of some parameters associated 
with them. This observable is sensitive both to the volume of space 
and to the growth of structure on scales of the order of a few tens of 
Mpc, as a function of redshift, more specifically, to the expansion 
history of the Universe and the amplitude of the primordial density 
perturbations on those scales. However, even slight deviations from 
Gaussianity in the probability distribution of the primordial density 
perturbations can have a measurable impact on the galaxy cluster 
abundance, especially at hight redshift and masses (see |Matarrese| 
|et al . 2000 Sef usatti et al.|200Tl and references therein). 

The number of galaxy clusters in bins of redshift, centered at 
redshift z, and with width Az, can be computed as 



N{z,±Kz) = Uy 
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where /j,!.,, is the fraction of sky being observed, dn/dM is the halo 
mass function and dV/dz is the comoving volume element, which 
in the case of a flat cosmology, is given by 
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with;^' (z) being the comoving radial distance. 

We will use the same prescription as in [Trindade et al.| 
( [201 2) to calculate the halo mass function in the presence of non- 
Gaussisuiities in the probability distribution of the primordial den- 
sity perturbations, which in turn followed[Lo Verde et al.|((2008[). 
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with a = 0.707, A = 0.322184 and p 
mass function ([Sheth et al.|2001jl, and 
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where S^m is the skewness of the smoothed density field [if /nl = 
0, then SsM = and Eq. (pj reduces to the Gaussian halo mass 
function]. This later formula comes from an extension to non- 
Gaussian density fields of the Press - Schechter formalism ( [Press[ 
[&Schechter[1974[ >, proposed by [Lo Verde et al.[ l [2008") . It has been 



shown to provide a good fit to results from N-body simulations (see 
[Wagner et al.|2010[ and references therein). 

The present-day value of the variance of the density perturba- 
tions, Sr, smoothed on a scale R can be computed as. 



, f d^ki r d^k2 
J (2nf J (iKf 
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(6) 



while the three-point function for the smoothed density field is 
given by ( |Lo Verde et al.|2008^ 

with (i = f(kr), Ti = W{ki,R}M( ki)T{ki), where T jk ) 
is the transfer function adopted from Bardeen et al. (1986|l, 
M{k) = (2/5) c^n;„^Hg-k^, W(k,R) is the smoothing top-hat win- 
dow and we use the shape parameter given by [Sugiyama| ( |1995[ ), 
r = n,„/!exp[-Qfc (l + V2A/fi„)]. 



3 RESULTS 

In order to quantify the biases which may arise in the estimation of 
cosmological parameters using the redshift evolution of the galaxy 
cluster abundance due to wrongly assuming the absence of primor- 
dial non-gaussianities, we will follow the same approach consid- 
ered in Sahlen e t al.[ ( |2009l >. Therefore, we will assume that the 
likelihood, X, of observing a given number of clusters, for a certain 
combination of cosmological parameters, in each bin of redshift 
and mass is given by 



InX: 



«; I'M 

ZZI 



A^^lnA'^- 



■ N'. 



inr^ + i)] 



(8) 



where n- and n^ are respectively the number of redshift and mass 
bins, while A'/, and N'^" are the number of counts for respectively 
the fiducial combination of cosmological parameters and the com- 
bination assumed to be true, for the i-th redshift bin and the j-th 
mass bin. Thus, our assumed observed cluster catalogue will not 
be a particular realization of the cluster redshift distribution in the 
context of the fiducial model, but a sort of average-catalogue. Nev- 
ertheless, this will allow for a very good estimate of the size and 
shape of the expected likelihood contours expected from an anal- 
ysis of the information contained in a real cluster catalogue, and 
avoid the offset in the best fit away from the fiducial values for the 
cosmological parameters, that would result from using a randomly- 
generated cluster catalogue based on the fiducial model (for more 
details see [Sahlen et al.|2009| l. 

In the above expression for the likelihood, we assumed that all 
clusters are randomly distributed in space, i.e. their positions follow 
a Poisson distribution. However, those positions are in fact spatially 
correlated. The effect of these correlations on the cluster abundance 
is often referred to as cosmic variance, while the variations in the 
cluster abundance due to the cluster distribution being a Poisson 
process are known as shot-noise. It has been shown, assuming pri- 
mordial Gaussian density perturbations, that the contribution of the 
cosmic variance to the statistical uncertainty associated with the 
galaxy cluster abundance is at least an order of magnitude smaller 
than the contribution due to shot-noise, almost independently of the 
surveyed sky area, as long as the cluster mass threshold is above 
5 X 10'**/!"' Mo IValageas et al.||201l| . Therefore, so that we can 
safely neglect the contribution from cosmic variance, we will set 
the cluster mass threshold to 5 x IO'^/i^'Mq, noting that the level 
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Figure 1. The most probable values of Clje, w and erg, and associated Icr (inner thin ticks) and Icr (outer thick ticks) confidence levels, for a sky coverage of 
4000 deg2. 



of non-Gaussianity we are assuming is relatively small, not affect- 
ing much the galaxy cluster abundance with respect to the Gaussian 
case jTrindade et al.|20T2) >. 

We assume a flat ACDM fiducial cosmology (i.e. the fidu- 
cial value for w is assumed to be -1), as derived in [Planck Col-| 
[laboration et al. | pOTSal, namely, a Hubble constant, Hq, equal 
to 100/ikms"' Mpc" with h = 0.673, fractional energy densities 
associated with dark energy and baryons today of Q^^ = 0.685 
and Q.ijh^ = 0.02205 respectively, a scalar spectral index, n,, equal 
to 0.9603, and power spectrum normalization of erg = 0.829. 
The level of non-gaussianity is parametrized by the /^.^ parameter, 
which we allow to vary from - 100 up to 100 with increments of 50. 
We also consider a fiducial sky area of 4000 deg^ and we generate 
the expected number of clusters in redshift bins with width Az = 0. 1 
up to a redshift of 1 .4. The cosmological parameters whose values 
we attempted to constrain using the information contained on those 
catalogues were O^^, w and erg, where w is the (assumed indepen- 
dent of redshift) constant equation of state associated with the dark 
energy. Flat priors were associated with all, rendering their pos- 
terior probabilities proportional to the likelihood given by Eq. [8] 
We set /jvL = as a prior, given our objective of quantifying the 
biases which may arise due to wrongly assuming the absence of 
primordial non-gaussianities. The exploration of the likelihood in 
the defined parameter space was carried out using a custom code 
based on standard Monte Carlo Markov Chain techniques (e.g. |von| 
|Toussaint|2011| ). 

As can be seen in Fig. 1, the biases on Clje and w, that arise 
from wrongly assuming f^t to be zero, are very small. In fact, it 
would be necessary larger values of 1/^.^1 than the ones considered 
here, combined with multiple mass bins and a larger sky coverage 
for the fiducial values of Cli/g and/or w to fall outside the derived 2cr 
confidence levels. 

Contrary to fi^^ and w, the bias on the estimation of erg due 
to wrongly assuming gaussianity is significantly more severe, al- 
though even in this case the exclusion of the fiducial value at more 
than 2cr requires a value for ff^L lower than about -80, or signif- 
icantly higher than 100. The same exclusion level is attained for 
lower l/jvil if multiple mass bins are considered and/or the sky area 
coverage increased. In the later case, we have confirmed that as ex- 
pected, the uncertainty associated with the estimation of each pa- 
rameter is inversely proportional to the square root of the sky area 
coverage. For example, increasing it to 40000 deg^ would result in 
the fiducial value for ag assumed here to be excluded at more than 
2cr for /^-i smaller than -30. 



The best fit values of Q.^,,, w and erg as a function of the non- 
Gaussianity level are given approximately by 



w 
erg 



0.686 + 3.3 X \0-''fNL + 3.7 x 10"^ {fNtf (9) 

-1.010 + 8.9 X IO^Va-l - 7.4 x lO"** (/wz,)- (10) 

0.829 + 1.7 X 10"V~L - 3.2 x 10"^ {fNtf , (H) 



with a maximum fitting error below 1%. The induced systematic er- 
rors in the cosmological parameters due to a systematic error in the 
non-Gaussianity parameter can be obtained by simply differentiat- 
ing Eqs. 9 to 1 1 with respect to /nl- Although the results above are 
specifically for the local parametrization of /jvl, they can be used 
to infer the what would happen had we considered the equilateral 
parametrization, given that a value for f^''^' "'"" which is 3 times 
some Z^™"' will yield roughly the same cluster abundance. We have 
also found that the dependences expressed through Eqs. 9 to 1 1 do 
not change significantly if the cluster mass threshold, assumed here 
to be 5 X 10'^/i^'Mq, is changed. 

We have found that the abundance of galaxy clusters can also 
be computed for the equilateral parametrization using 



rlocal 

Jnl 



requil 

Jnl 
3.6 



(12) 



which seems to be consistent with the recent results of lShandera et| 
|ar|j2013]> and do not change much with cluster mass threshold. 



4 CONCLUSIONS 

In this paper we have computed the dependence, as a function of 
/ail, of the biases that arise in the estimation of the cosmological 
parameters Q^c, w and erg, when it is (wrongly) assumed, for the 
purpose of the statistical analysis, that the density field is Gaussian 
distributed. We have found that such biases are quite small for the 
first two parameters, but significant in the case of erg, in particular 
in the face of the high statistical accuracy with which this parameter 
is expected to be determined in the near-future (e.g. |Pillepich et af] 
[20T2l >. 

If /ivL is assumed to be scale-independent, then the recent 
results obtained by the Planck team ( [Planck Collaboration et| 
[al. |[2013F^ , impose severe constraints on the amount of non- 
Gaussianity at cluster scales {f^L = 2.7 ± 5.8 for the local param- 
eterization) making it safe to neglect primordial non-Gaussianities 
in the determination of cosmological parameters using the galaxy 
cluster abundance. However, fjyi could be a function of scale (see 
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|Lo Verde et al.|2008l and references therein) , in which case the ef- 
fects of non-gaussianity on the cluster abundance may need to be 
taken into account in the determination of erg . In fact, Eq. (11) sug- 
gests that a significantly negative value for f^t at cluster scales, of 
the order of -240, would reconcile the constraint on erg obtained 
using the cluster abundance inferred from the Planck Sunyaev- 
Zeldovich Catalogue {a^^ = 0.11 ± 0.02, [Planck Collaboration et| 
|al. |2013c) with the one obtained from the properties of the primor- 
dial CMB temperature anisotropies (erg = 0.829 ± 0.012, [Pianck| 
ICollaboration et al. |20l3al i. 



ACKNOWLEDGEMENTS 

We thank Lara Sousa, Antonio da Silva and Michael Bazot 
for useful discussions during the preparation of this paper. 
AMMT was supported by the FCT/IDPASC grant contract 
SFRH/BD/51647/2011. AMMT, PPA and PTPV acknowledge fi- 
nancial support from project PTDC/FIS/1 11725/2009, funded by 
Fundagao para a Ciencia e a Tecnologia. 



REFERENCES 



P. A. R., 



P. A. R., 



P A. R., 



Aghanim, N., 
Aghanim, N., 
Aghanim, N., 



et al. 2013a, 



et al. 2013b, 



et al. 2013c, 



Planck Collaboration, Ade, 

arXiv:1303.5076 
Planck Collaboration, Ade, 

arXiv:1303.5084 
Planck Collaboration, Ade, 

arXiv:1303.5080 
Pillepich, A., Porciani, C, Reiprich, T. H. 2012, MNRAS, 422, 44 
Shandera, S., Mantz, A., Rapetti, D., & Allen, S. W. 2013, 

arXiv:1304.1216 
Trindade, A. M. M., Avelino, P P, & Viana, P T. P 2012, MNRAS, 424, 

1442 
Hazra, D. K., & Guha Sarkar, T. 2012, Physical Review Letters, 109, 

121301 
Achitouv, I. E., & Corasaniti, P. S. 2012, J. Cosmology Astropart. Phys., 

2,2 
Tashiro, H., & Ho, S. 2012, arXiv: 1205.0563 

Takeuchi, Y, Ichiki, K., & Matsubara, T. 2012, Phys. Rev. D, 85, 043518 
Maio, U., Salvaterra, R., Moscardini, L., & Ciardi, B. 2012, MNRAS, 426, 

2078 
Giannantonio, T, Porciani, C, Carron, J., Amara, A., & Pillepich, A. 

2012, MNRAS, 422, 2854 
Schafer. B. M., Grassi, A., Gerstenlauer, M., & Byrnes, C. T. 2012, MN- 
RAS, 421, 797 
Hilbert, S., Marian, L., Smith, R. E., & Desjacques, V. 2012, MNRAS, 

426, 2870 
Alishahiha, M., Silverstein, E., & Tong, D. 2004, Phys. Rev. D, 70, 123505 
Avelino, P P & Viana, P T. P 2000, MNRAS, 314, 354 
Baldauf, T, Seljak, U., Senatore, L., & ZaldaiTiaga, M. 2011, ArXiv e- 

prints 
Bardeen, J. M., Bond, J. R., Kaiser, N., & Szalay, A. S. 1986, ApJ, 304, 

15 
Wang, S., Khoury, J., Haiman, Z., & May, M. 2004, Phys. Rev D, 70, 

123008 
Chiu, W. A., Ostriker, J. P, & Strauss, M. A. 1998, ApJ, 494, 479 
Creminelli, F 2003, J. Cosmology Astropart. Phys., 10, 3 
Fedeli, C, Moscardini, L., & Matan-ese, S. 2009, MNRAS, 397, 1125 
Fosalba, P., Gaztanaga, E., & Elizalde, E. 2000, in Astronomical Society 

of the Pacific Conference Series, Vol. 200, Clustering at High Redshift, 

ed. A. Mazure, O. Le Fevre, & V. Le Brun. 408-+ 
Komatsu, E. 2010, Classical and Quantum Gravity, 27, 124010 
Komatsu, E., Dunkley, J., Nolta, M. R., et al. 2009, ApJS, 180, 330 
Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18 
Komatsu, E. & Spergel, D. N. 2001, Phys. Rev D, 63, 063002 



Lo Verde, M., Miller, A., Shandera. S., & Verde, L. 2008. J. Cosmology 

Astropart. Phys., 4, 14 
Lyth, D. H. & Rodriguez, Y. 2005, Physical Review Letters, 95, 121302 
Maldacena, J. 2003, Journal of High Energy Physics, 5, 13 
Matarrese, S. & Verde, L. 2008, ApJ, 677, L77 
Matarrese, S., Verde, L., & Jimenez, R. 2000, ApJ, 541, 10 
Press, W. H. & Schechter, P 1974, ApJ, 187, 425 
Riotto, A. & Sloth, M. S. 2011, Phys. Rev D, 83, 04I301(R) 
Robinson, J. & Baker, J. E. 2000, MNRAS, 311, 781 
Robinson, J., Gawiser, E., & Silk, J. 2000, ApJ, 532, 1 
Sahlen, M., Viana, P T. P, Liddle, A. R., et al. 2009, MNRAS, 397, 577 
Salopek, D. S. & Bond, J. R. 1990, Phys. Rev. D, 42, 3936 
Seery, D. & Lidsey, J. E. 2005. J. Cosmology Astropart. Phys., 6, 3 
Sefusatti, E., Vale, C, Kadota, K., & Frieman, J. 2007, ApJ, 658, 669 
Sefusatti, E. & Komatsu, E. 2007, Phys. Rev D, 76, 083004 
Sheth, R. K., Mo, H. J., & Tormen, G. 2001, MNRAS, 323, 1 
Slosar, A., Hirata, C, Seljak, U., Ho, S., & Padmanabhan, N. 2008, J. 

Cosmology Astropart. Phys., 8, 31 
Sugiyama, N. 1995, ApJS, 100, 281 

Valageas, A., Clerc, C, Pacaud, U., & Pierre, N. 201 1, A&A, 536, A95 
von Toussaint, U. 201 1, Reviews of Modern Physics, 83, 943 
Wagner, C, Verde, L., & Boubekeur, L. 2010, J. Cosmology Astropart. 

Phys., 10, 22 



© 2013 RAS, MNRAS 000, ??-?? 



